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EVALUATION  OF  MARSH/ ESTUARINE  WATER  QUALITY  AND 
ECOLOGICAL  MODELS:  AN  INTERIM  GUIDE 

Purpose  and  Scope 

1.  This  report  introduces  marsh/estuarine  water  quality  and 

ecological  models  and  assists  in  evaluating  the  applicability  of  matne-  | 

matical  modeling  methodologies  to  environmental  assessments  in  estuaries 

and  intertidal  wetlands.  This  report  provides  preliminary  guidance  to 
District  planners,  engineers,  and  resource  managers. 

2.  A  brief  summary  of  water  quality  and  ecological  mod 'ling  in  ; 

the  marsh/estuarine  environment  is  provided  as  the  basis  for  the  model¬ 
ing  recommendations. 

Mater  Quality  and  Ecological  Problems  1 

3.  A  survey  of  Corps  of  Engineers  (CE)  field  offices  with  coastal  | 

responsibilities  identified  existing  and  anticipated  water  quality  and  \ 

ecological  problems  associated  with  CE  activities  in  coastal  marshes  and  1 

estuaries.  Existing  and  potential  problems  were  classified  and  tabulated  * 

as  to  marine  wetlands,  coastal  zone,  an^  environmental  monitoring  (Ta-  I 

ble  1).  The  feasibility  of  numerical  modeling  techniques  for  addressing  ; 

\ 

each  identified  problem  is  noted.  A  more  detailed  description  of  appli-  1 

cable  modeling  approaches  is  outlined  in  Table  2.  5 

Marine  wetlands 

A.  Mathematical  modeling  methodologies  are  not  sufficiently  devel¬ 
oped  to  address  the  marine  wetlands  problems.  Existing  research  io 
quantifying  material  cycling  and  primary  and  secondary  production  within 
the  marsh  and  between  the  marsh  and  adjacent  coastal  waters.  Existing 
mathematical  models  addressing  these  phenomena  are  research  tools  and 
are  not  suitable  for  planning  and  management. 

5.  Conceptual  models  should  be  used  to  guide  studies  that  provide 
increased  understanding  and  quantification.  The  potential  applicability 
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of  marsh  and  estuarine  ecosystem  models  to  these  problems  suggests  that 
studies  be  designed  to  understand  important  processes  and  to  provide 
data  for  future  model  verification. 

Coastal  zone 

6.  Mathematical  modeling  techniques  are  not  suitable  for  coastal 
zone  problems  12  and  13  (Table  1).  However,  mathematical  models  can  be 
applied  to  problems  14-16  when  used  in  conjunction  with  complementary 
study  procedures. 

7.  Present  mathematical  modeling  studies  address  the  fate  of 
dredged  material.  Mathematical  modeling  to  assess  the  environmental 
effects  of  turbidity  resulting  from  dredging  and  disposal  is  not  recom¬ 
mended  at  this  time  (problem  14). 

8.  Problems  15  and  16  are  directly  addressable  through  mathe¬ 
matical  modeling.  However,  a  more  detailed  understanding  of  the  hydro¬ 
dynamics  of  very  shallow  waters — including  sheet  flow  in  marshes;  flow 
in  very  shallow  channels;  sediment  transport;  periodic  inundation;  and 
the  relationships  of  residence  time  to  winds,  tides,  and  freshwater  in¬ 
flow  and  rainfall — is  needed.  Present  understanding  of  biological  and 
chemical  processes  would  limit  simulations  to  assessing  excessive  algal 
growth  and  dissolved  oxygen  relationships.  Additional  knowledge  would 
be  required  t.o  simulate  the  cycling  of  toxic  constituents. 

Environmental  monitoring 

9.  Mathematical  modeling  methodologies  are  not  appropriate  to 
address  environmental  monitoring  problems. 

Marsh/ Estuarine  Water  Quality  and  Ecological  Models 

10.  District  planners,  engineers,  and  resource  managers  require  an 
understanding  of  the  behavior  of  marsh/estuarine  ecosystems  and  how  the 
system  would  respond  as  a  result  of  future  CE  activities. 

11.  One  way  to  assess  the  effects  of  proposed  activities  is  to 
make  the  actual  changes  in  the  system  and  directly  observe  the  effects; 
however,  such  direct  experimentation  is  usually  impractical.  An  alter¬ 
native  is  to  model  the  physical  system  in  such  a  way  that  essential 


cause-and-effect  relationships  are  maintained.  The  model  may  be  a 
scaled  physical  representation  or  a  mathematical  description  based  on 
physical  laws  and  empirical  formulae. 

12.  In  a  mathematical  model,  the  cause-and-effect  relationships 
between  system  elements  are  represented  by  a  set  of  mathematical  equa¬ 
tions.  Mathematical  models  may  vary  from  a  few  equations  that  can  be 
solved  by  hand  computation  to  hundreds  of  equations  requiring  the  use  of 
a  digital  computer.  The  complexity  of  the  mathematical  model  selected 
depends  largely  on  the  questions  to  be  asked  about  the  behavior  of  the 
actual  system. 

13.  Marshes  and  estuaries  are  complex  Interactive  ecosystems. 
Salient  features  that  contribute  to  their  complexity  are  diverse  topog¬ 
raphy,  environmental  gradients,  and  the  large  number  of  biological 
components. 

1A,  No  single  model  exists  that  can  simulate  all  aspects  of  water 
quality  and  ecological  processes  in  marsh/estuarine  ecosystems.  If  such 
a  model  were  developed,  the  resulting  complexity,  extensive  data  require¬ 
ments,  and  cost  of  computation  would  make  its  use  impractical.  An 
alternative  strategy  is  a  wide  variety  of  more  specialized  models  devel¬ 
oped  to  more  effectively  address  particular  aspects  of  water  quality  or 
ecosystem  processes.  This  section  of  the  report  will  review  the  types 
and  uses  of  mathematical  models. 

Hydrodynamic  models 

15.  The  basis  for  water  quality  models  and  ecosystem  models  that 
use  the  principle  of  conservation  of  mass  is  a  model  of  the  circulation 
of  the  estuary.  The  circulation  in  a  partially  mixed  estuary  is  mainly 
driven  by  the  tide,  river  inflow,  salinity  gradient,  and  wind.  It  is 
now  considered  that,  in  many  estuaries,  all  four  forcing  mechanisms  are 
equally  important  in  influencing  the  circulation  and  the  transport  of 
water  quality  constituents.  Estuarine  circulations  are  some  of  the  most 
complex  geophysical  flows  in  existence. 

16.  Physical  oceanography  of  estuaries  is  discussed  in  the  texts 
by  Dyer  (1973),  Officer  (1976),  and  McDowell  and  O'Connor  (1977);  sym¬ 
posia  edited  by  Kjerfve  (1978)  and  Hamilton  and  Macdonald  (1980);  and 
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reviews  by  Liu  and  Leendertse  (1978)  and  Hamilton  and  Fucik  (1980). 

17.  One-dimensional  models.  One  use  of  one-dimensional  models 
that  included  an  advection-dispersion  equation  was  the  study  of  salinity 
intrusion  in  an  estuary  (Stigter  and  Siemons  1967;  Thatcher  and  Harlaman 
1972;  Williams  and  West  1973).  Neglect  of  the  variation  with  depth  of 
current  and  salinity  meant  that  the  upstream  flux  of  salt  due  to  density 
currents  and  mixing  due  to  vertical  phase  differences  in  the  tidal 
variations  of  current  and  salinity  were  lumped  together  in  a  one¬ 
dimensional  longitudinal  dispersion  coefficient.  Because  the  dispersion 
coefficient  varied  among  estuaries  and  was  highly  dependent  upon  river- 
flow  and  distance  from  the  estuary  mouth  (Williams  and  West  1973),  one¬ 
dimensional  models  were  difficult  to  use  in  a  predictive  mode  in  par¬ 
tially  mixed  estuaries.  However,  these  models  serve  as  the  hydrodynamic 
basis  for  the  majority  of  estuarine  water  quality  models. 

18.  Two-dimensional  models.  Consideration  of  the  depth  dimension 
permits  the  inclusion  of  the  majority  of  the  fundamental  physical  pro¬ 
cesses  that  affect  the  advection  and  dispersion  of  water  quality  con¬ 
stituents.  The  first  study  demonstrating  the  importance  of  depth  varia¬ 
tion  on  currents  and  salinities  was  the  steady-state  analytical  solution 
of  Rattray  and  Hansen  (1962) .  These  equations  were  expanded  to  explic¬ 
itly  include  the  longitudinal  salt  balance  (Hansen  and  Rattray  1965; 
Rattray  1967;  Winter  1973;  Festa  and  Hansen  1976). 

19.  Time-dependent  equations  for  a  partially  mixed  estuary  were 
numerically  solved  by  Hamilton  (1975)  and  Blumberg  (1975).  The  solution 
techniques  were  explicit  with  time  steps  of  about  1  min.  More  recent 
semi- implicit  methods  permit  time  steps  of  about  30  min  (Hamilton  1976; 
Wang  and  Kravitz  1980),  thereby  permitting  wind-driven  transients  to  be 
investigated  (Wang  1980) . 

20.  The  vertically  averaged  equations  of  motion  have  been  exten¬ 
sively  used  for  storm  surge  modeling.  Hinwood  and  Wallis  (1975a,  1975b) 
provide  a  review  of  these  models.  Vertically  averaged  hydrodynamic 
models  are  now  quite  sophisticated  employing  variable  grids  and  embedded 
subgrid  features  to  handle  complex  geometry  such  as  channels  and  barriers 
(Butler  1980;  Reid  et  al.  1977;  Reid,  Vastano,  and  Reid  1977). 
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21.  A  few  investigators  have  included  vertically  integrated 
advection-dispersion  equations  for  salinity  and  other  water  quality 
constituents.  Leendertse  and  Gritton  (1971)  hindcasted  water  quality 
constituents  in  Jamaica  Bay,  N.  Y.,  and  Hess  and  White  (1977)  modeled 
the  dispersal  of  a  marked  fluid  released  in  Narragansett  Bay,  R.  I._ 

22.  Three-dimensional  models.  Three-dimensional  models  are  very 
complex  and  only  a  few  have  been  developed.  The  model  by  Leendertse  and 
Liu  (1975,  1977)  uses  explicit  finite  difference  methods  and  thus  is 
extremely  costly  in  computer  time  for  any  lengthly  simulation.  Capon! 
(1976)  developed  a  time-dependent  model  that  used  the  complete  three- 
dimensional  Navier-Stokes  equations  instead  of  the  Boissinesq  approxima¬ 
tion.  However,  in  an  application  to  Chesapeake  Bay,  the  use  of  a  rela¬ 
tively  coarse  grid  may  have  resulted  in  some  anomalous  circulations. 

Wang  (1980)  proposed  the  use  of  semi-implicit  methods  and  mode  splitting 
for  increasing  the  computational  efficiency  of  three-dimensional  models. 

23.  Summary.  Table  3  provides  examples  of  the  hydrodynamic  models 
discussed.  The  list  is  not  exhaustive,  but  is  intended  to  show  the 
variety  of  approaches  to  the  numerical  modeling  of  estuarine 
hydrodynamics . 

24.  Many  one-  and  two-dimensional  models  include  salinity  advec- 
tion  and  dispersion.  Some  models  use  tidally  averaged  velocity  as  input 
to  salinity  transport  models.  Investigations  by  Harleman  (1971)  indi¬ 
cate  differences  between  the  results  of  tidally  averaged  and  time- 
dependent  one-dimensional  simulations  of  salinity  transport.  Further¬ 
more,  he  concluded  that  the  importance  of  the  highly  empirical 
longitudinal  dispersion  coefficient  decreases  as  the  accuracy  of  descrip¬ 
tion  of  the  advective  motion  increases. 

25.  The  alternatives  to  using  tidally  averaged  or  net  tidal  veloc¬ 
ities  are  to  output  the  hydrodynamic  computational  results  for  subsequent 
input  to  salinity  transport  model  or  to  embed  the  salinity  transport 
equations  within  the  hydrodynamic  computations.  Embedding  permits  cou¬ 
pling  of  the  salinity  gradient  to  the  momentum  equations.  However, 
highly  discretized  hydrodynamic  models  have  inherent  stability  problems 
requiring  time  steps  on  the  order  of  seconds  to  minutes  with  maximum 
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simulation  periods  of  days,  while  water  quality  phenomena  of  interest 
such  as  salinity  distribution  and  phytoplankton  growth  may  require 
simulation  periods  on  the  order  of  weeks  or  months,  and  possibly  longer. 
Water  quality  and  ecological  models 

26.  Water  quality  models.  Historically,  the  analysis  of  water 
quality  has  concentrated  on  the  DO  and  BOD  due  to  waste  loads.  The 
majority  of  DO-BOD  models  evolved  from  the  study  by  Streeter  and  Phelps 
(1925),  who  assumed  that  the  balance  between  DO  and  BOD  concentrations 
was  the  result  of  two  processes:  the  reaeration  of  the  water  column  and 
the  consumption  of  DO  in  the  oxidation  of  BOD. 

27.  Later  modeling  emphasis  has  been  on  extending  and  refining  the 
Streeter-Phelps  formulation  by  using  a  more  generalized  mass  balance 
approach  and  by  the  inclusion  of  additional  processes  such  as  benthic 
oxygen  demand,  benthic  scour  and  deposition,  photosynthesis  and 
respiration  of  aquatic  plants,  and  nitrification  (Dobbins  1964;  O'Connor 
1967;  Dresnack  and  Dobbins  1968;  O'Connor  and  DiToro  1970;  Grantham, 
Schaake,  and  Pyatt  1971;  Whitehead  and  Young  1975). 

28.  The  more  comprehensive  estuarine  water  quality  models  have 
been  developed  to  include  the  nitrogen  and  phosphorus  cycle  and  the 
lower  trophic  levels  of  phytoplankton  and  zooplankton  (Chen  and  Orlob 
1972;  DiToro  et  al.  1977;  Dahl-Madsen  1978). 

29.  Water  quality  modeling  is  discussed  in  the  texts  edited  by 
lames  (1978);  Canale  (1976);  Biswas  (1976);  and  Hatzinger,  Vanlelgveld, 
and  Zoeteman  (1978);  and  in  review  articles  by  DiToro  et  al.  (1977); 

Hahn  and  Schreiner  (1978);  and  O'Connor,  Thomann,  and  DiToro  (1977). 

30.  A  few  water  quality  modeling  investigations  of  estuaries  have 
addressed  water  quality  constituents  other  than  DO-BOD,  nitrogen,  phos¬ 
phorus,  and  the  lower  trophic  levels  of  the  biological  community.  A 
number  of  investigators  have  modeled  the  algal  nutrient  silica  (Festa 
and  Hansen  1976;  Peterson,  Festa,  and  Conomos  1978;  Rattray  and 
Officer  1979). 

31.  Nihoul  et  al.  (1979)  and  Billen  and  Smitz  (1977)  modeled 
selected  chemical  constituents  in  an  estuary  by  assuming  thermodynamic 
equilibrium.  The  variable  controlling  partitioning  between  constituents 


8 


■V* ,  fi-^‘  -*v\tar  ««•  *&■**#  ^  S&^PR 

v  J’  *  -$ 


was  redox  potentials  estimated  by  considering  the  oxygen  budget.  The 
authors  admit  that  thermodynamic  equilibrium  is  not  necessarily  fulfilled 
throughout  the  estuary. 

32.  There  have  been  very  few  models  simulating  fate  and  effects  of 
heavy  metals.  The  fate  of  metals  in  estuaries  is  very  complicated  and 
unclear,  involving  adsorption-desorption  reactions,  flocculation,  precipi¬ 
tation,  sedimentation,  and  biological  uptake,  Jorgensen  (1979)  modeled 
heavy  metals  using  the  concept  of  trophic  length  lev^l.  The  concept 
involved  the  consideration  of  a  length  scale  introduced  into  the  equa¬ 
tions  as  an  independent  variable  that  represents  the  time  a  constituent 
remains  in  a  particular  trophic  level. 

33.  Ecological  models.  In  contrast  to  water  quality  models,  eco¬ 
logical  models  are  more  descriptive,  emphasizing  exhaustiveness  and 
resolution  often  to  the  limits  of  potential  data  availability.  Ecologi¬ 
cal  models  include  numerous  biological  species  or  species  aggregates  and 
emphasize  food  chain  and  species  interactions. 

34.  Cushing  (1975),  in  his  book  on  marine  ecology  and  fisheries, 
points  out  that,  in  marine  ecological  studies,  the  limiting  nutrient  con¬ 
cept  has  been  abandoned  with  an  increased  emphasis  placed  on  predator- 
prey  relationships.  Early  marine  ecosystem  models  concentrated  on  the 
nutrient-phytoplankton-zoon1 ankton  interactions  with  complex  formula¬ 
tions  describing  the  transfer  of  material  between  the  components  (Steele 
1974).  Early  investigations  did  not  include  explicit  spatial  variations, 
but  more  recent  studies  allow  loth  spatial  and  temporal  variability  in 
systems  which  have  large  physical  gradients  such  as  upwelling  regions 
(Walsh  1971;  O’Brien  and  Wroblewski  1972)  and  estuaries  (Winter,  Banse, 
and  Anderson  1975;  Kremer  and  Nixon  1978). 

35.  In  contrast  to  the  early  marine  models  with  few  components  but 
complex  interactions  is  the  approach  advocated  by  Patten  (1968). 

Patten's  approach  uses  very  complex  model  structures  that  involve  many 
species,  but  the  interactions  between  components  are  generally  simple, 
linear,  and  empirical.  These  models  tend  to  be  very  theoretical  since 
little  data  are  available  for  calibration. 

36.  Estuarine  ecological  models.  Few  estuarine  erosytem  models 
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exist,  Pomeroy  et  al.  (1972)  simulated  phosphorus  flux  in  several 
Georgia  estuaries.  Bahr  (1974)  simulated  energy  flux  in  Georgia  oyster 
beds.  Dame,  Vernbert,  and  Bonnell  ( J 9 7 7 )  used  a  similar  model  for  South 
Carolina  oyster  beds  emphasi2ing  sensitivity  analysis.  Show  (1979) 
developed  a  hydrological-biological  model  to  investigate  zooplankton 
distributions  in  a  Texas  coast  embayment.  Ferguson  and  Adams  (197?) 
simulated  epifauna  and  juvenile  fish  of  an  eelgrass  community  to  investi¬ 
gate  their  response  to  temperature. 

37.  The  two  most  comprehensive  ecological  models  of  estuaries  are 
those  by  Winter,  Banse,  and  Anderson  (1975)  and  Kremer  and  Nixon  (1978). 
Winter,  Banse,  and  Anderson  simulated  primary  production  for  the  central 
basin  of  Puget  Sound,  which  is  a  deep  fjord  with  strong  stratification. 
The  model  of  Kremer  and  Nixon  was  applied  to  Narragansett  Bay,  which  is 
a  shallow,  wide,  and  partially  mixed  estuary.  Both  models  were  based  on 
extensive  investigations  by  many  researchers  over  many  years, 

38.  Marsh  ecological  models.  Very  few  marsh  ecological  models 
exist.  Wiegert  and  Wetzel  (1979)  developed  a  model  of  a  Georgia  Spartina 
marsh.  The  objective  of  model  development  was  to  describe  the  pathways 
and  dynamics  of  carbon  flow  in  the  marsh.  Hopkinson  and  Day  (1977)  simu¬ 
lated  carbon  and  nitrogen  flows  in  a  Louisiana  salt  marsh. 

39.  Models  have  been  developed  to  address  specific  components  of 
the  marsh.  Nixon  and  Oviatt  (1973)  simulated  diurnal  DO  in  tidal  creeks 
and  embayments  to  investigate  the  effects  of  sewage  input  and  temperature 
increases  on  the  DO  concentrations  in  the  embayments.  Reimold  (1974) 
specifically  addressed  the  effect  on  the  Spartina  system  of  perturbations 
to  the  marsh.  Ziernan  and  Odum  (1978)  simulated  plant  growth  and  succes¬ 
sion  on  an  estuarine  salt  marsh.  Lugo,  Sell,  and  Snedaker  (1976)  simu¬ 
lated  mangrove  production. 

40.  Summary .  Table  4  provides  a  summary  of  the  different 
approaches  to  water  quality  modeling.  The  list  is  not  exhaustive  but 
provides  examples  ranging  from  one-dimensional,  tidally  averaged  DO-BOD 
models  through  two-dimensional,  time-dependent  models  that  include  DO- 
BOD,  plant  nutrients,  and  the  lower  trophic  levels. 

41.  Table  5  summarizes  water  quality  models  modified  or  developed 
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through  grants  by  State  and  Federal  water  resource  agencies  and  are  thus 
readily  available  with  user  documentation. 

42.  The  SEM  and  ES001  models  provide  the  one-dimensional,  tidally 
averaged  analytical  solution  initially  discussed  by  O'Connor  (1960). 

The  DEM  model,  and  the  derived  models  RECEIV  and  RIVSCI,  are  based  on 
the  WRE  Link-Node  hydrodynamic  model.  The  model  is  basically  a  one- 
dimensional  solution  adapted  to  two-dimensional  estuaries.  The  model 
has  been  characterized  as  a  good  descriptive  tool  but  lacks  a  predic¬ 
tive  capability.  The  ESTECO  model  was  developed  for  the  Texas  Water 
Development  Board  and  uses  a  two-dimensional,  vertically  averaged  hydro- 
dynamic  model  called  HYDTID  to  generate  a  velocity  field  for  the  water 
quality  module.  The  water  quality  module  is  based  upon  the  concepts 
originally  included  in  the  DEM  series  of  models. 

43.  In  summarizing  the  review  of  marsh/estuarine  water  quality 
and  ecological  models,  the  following  conclusions  can  be  drawn: 

a.  Additiorai  studies  are  needed  to  investigate  the  appro¬ 
priate  procedures  to  couple  hydrodynamic  and  water 
quality  models  since  the  time  scales  of  interest  may 
differ. 

b.  The  DO-BOD  estuarine  models  are  appropriate  for  investi¬ 
gating  the  assimilative  capacity  of  estuaries  to  heavy 
organic  enrichment.  However,  for  mildly  perturbed  estu¬ 
aries  in  which  photosynthesis,  algal  respiration,  decom¬ 
position,  and  mixing  processes  play  dominant  roles,  the 
understanding  and  characterization  of  significant  pro¬ 
cesses  are  less  well  known. 

c.  Models  developed  to  address  eutrophication  due  to  nutrient 
enrichment  have  not  been  verified  in  most  cases.  These 
models  are  not  capable  of  predicting  absolute  values  under 
varying  environmental  conditions.  These  models  may  be  use¬ 
ful  in  evaluating  minor  perturbations,  such  as  increased 
nutrient  loadings  or  turbidity.  However,  these  perturba¬ 
tions  cannot  be  catastrophic.  Events  that  significantly 
alter  the  species  composition  presently  cannot  be  eval¬ 
uated  through  simulations. 

d_.  Detailed  mathematical  ecological  models  presently  are  not 
sufficiently  developed  to  be  applicable  to  environmental 
impact  analysis.  The  extensive  data  requirements  exceed 
the  capabilities  of  most  specific  project  studies. 
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Table  1 

Existing  and  Potential  Water  Quality  and  Ecological  Problems 
Identified  by  a  Survey  of  CE  Field  Offices  with 
Coastal  Responsibility 


Problem 

Description 

Modeling 

Feasibility* 

Marine  Wetlands 

1 

Existing  wetlands  need  to  be  surveyed. 

— 

2 

Problems  are  encountered  when  applying  the 
different  kinds  of  available  classifica¬ 
tion  schemes. 

" 

3 

Information  is  needed  on  the  role  of  buffer 
zones  around  wetlands. 

— 

A 

Wetland  values  need  to  be  assessed  by 
productivity  measurements. 

— 

5 

More  information  is  needed  on  the  manage¬ 
ment  of  diked  impoundments, 

— 

6 

More  information  is  needed  on  the  assimi¬ 
lation  of  effluent  and  dredged  material 
by  the  wetland. 

7 

Delta  growth  as  a  means  of  marsh  creation 
needs  to  be  investigated. 

— 

8 

Additional  understanding  of  marsh  creation 
through  dredged  material  disposal  is 
needed. 

9 

Importance  of  created  marsh  on  the  total 
lay-marsh  ecosystem  requires  additional 
study. 

10 

Effects  of  changes  in  freshwater  flow  on 
wetland  areas  need  to  be  investigated. 

— 

11 

Procedures  to  evaluate  the  effects  of 
stresses  on  the  wetlands  such  as  plat¬ 
forms  and  highway  bridges  need  to  be 
developed. 

(Continued) 

*  Entries  in  this  column  are  defined  as  follows:  —  mathematical 
modeling  methodologies  are  not  sufficiently  developed  to  address  the 
problem;  /  modeling  is  feasible. 
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Table  1  (Concluded) 

Problem 

Description 

Modeling 

Feasibility 

Coastal  Zone 

12 

Offshore  disposal  sites  are  needed  due  to  a 
lack  of  upland  sites. 

— 

13 

Problems  are  encountered  with  bioassays. 

— 

14 

Concerns  exist  over  the  effects  of  turbid¬ 
ity  produced  by  dredging  and  disposal 
operations. 

/ 

15 

Water  quality  problems  such  as  eutrophica¬ 
tion  or  pollution  exist  in  estuaries  or 
rivers. 

✓ 

16 

Procedures  for  evaluating  the  impacts  of 
cumulative  development  of  dead-end  canals, 
small  boat  harbors,  and  marinas  are 
needed. 

/ 

17 

Field  offices  need  hydrodynamic  and 
advection-dispersion  models. 

/ 

Environmental  Monitoring 

18 

Methods  of  storing  and  handling  data  and 
literature  for  environmental  assessments 
need  to  be  modernized. 

— 

19 

Longer  lead  times  are  needed  for  compre¬ 
hensive  studies  for  environmental  impact 
assessment. 

20 

More  long-term  monitoring  is  needed  to  per¬ 
mit  impact  evaluations. 

— 

21 

Long-term  studies  are  needed  to  provide 
additional  understanding  of  problems 
unique  to  the  field  office. 

22 

Exchange  of  technical  information  be¬ 
tween  field  offices  needs  to  be  improved. 

— 

Table  3 


Representative  Types  of  Hydrodynamic  Models  of  Estuaries* 


Dimensions 

Solution  Technique 

Variables 

Equations 

Forcing 

Application 

Elliott  (1975) 

X 

(2  Layer) 

Noncoupled  kinematic 
and  dynamic  box  model 

Vi<x) 

Uu,lW 

Continuity 

Salt  balance 

River flow 

Potomac  Estuary 

VS, 

u-m omen turn 
(approximated) 

Hamilton  (1975) 
Blumborg  (1975) 

Wang  and  Kravitz 
(1980) 

x,z,t 

Explicit 

Finite  differences 
Semi-implicit 

;(x,t) 

u(x,z,t) 

v(x,x,t) 

s(x.l,t) 

Continuity 
u -momentum 

Salt  conservation 

River flow 

A  tide 
+  wind 

Rotterdam  Waterway 

Potomac  Estuary 

Hodgins  (1979) 

x,t 

(2  layer) 

Finite  differences 
Scmi-impllcit 

;(*.t> 

h(x,t) 

uu,l<x>t) 

“u,l<X,t> 

Continuity 
u -momentum 

Salt  conservation 

Wind 

River flow 

Albernl  inlet, 
British  Columbia 

Johns  (1978) 

x,z,t 

Finite  differences 
Explicit  x-coord 
Semi-implicit 
z-coord 

;(x,t) 

u(x,i,t) 

v(x,Z,t) 

E(x,t,t) 

Continuity 
u-momentum 
Turbulent  energy 

Tide 

Rectangular 
uniform  channel 

Reid,  Vastano,  and 
Reid  (1977) 

x,y,t 

Explicit 

Finite  differences 
(subgrid  scale) 
(channels  &  barriers) 

C(x,y,t) 

u(x,y,t) 

v(x,y,t) 

Continuity 

u-momentum 

v-momentum 

Tide 

A  wind 

Sabine-Calcasieu 
area  of  gulf 
coast 

Pearson  and 

Winter  (1977) 

x,y,t 

Finite  elements 

(time  decomposed  Into 
tidal  harmomos  w) 

C(x.y.u) 

u(x,y,u) 

v(x,y,u) 

Continuity 

u-momentum 

v-momentum 

Tide 

Hood  Canal  Wash. 
(Jamart  personal 
communication) 

Liu  and 

Leendcrtse  (1978) 

x«y,z,t 

Finite  differences 
Semi-implicit 

t(x,y,t) 

u(x,y,z,t) 

v(x,y,z,t> 

<v(x,y,z,t) 

s(x,y,t,t) 

E(x,y,a,t) 

Continuity 

u-momentum 

vmomentun 

Salt  conservation 
Turbulent  energy 

Tide 

Wind 

Rlverflow 

Long  Island  Sound 
San  Francisco  Bay 

Chesapeake  Bay 

Butler  (1978) 

x.y.t 

Seml-lttpllclt 

(variable  grid,  sub¬ 
grid  scale  features) 

t(x.y.t) 

u(x,y,t) 

v(x,y,t) 

Continuity 
u -momentum 
v-momentum 

Tide 

Wind 

Rlverflow 

Tsunami 

Ma.onboro  Inlot  t 
estuary  (tide) 
Cor. on  Inlet  8 
eatuary  (tide) 

Great  Eg*  Harbor 
Inlet  A  estuary 
(tide) 

Coos  Bay  A  South 
Slough,  Oregon 
(tide) 

Bolsa  Chlca  Bay 
(tide) 

Los  Angeles  A 
Long  Beach 
Harbor  (tide) 
Galveston  Bay 
(tide  &  surge) 
Crescent  City,  Cal¬ 
ifornia  (tsunami) 
Hauula-Punsluu  re¬ 
gion,  Oahu,  Ha¬ 
waii  (tsunami) 


Variable 

x 

y 

Z  4 

t 

u 

u,v,w 

u,v 

s 

E 

A 

h 

VS, 

Superscripts  u,l 


_ Definition 

Longitudinal  horizontal  axis 
Across-channel  horizontal  axis 
Depth  axis 
Tine 

Harmonic  decomposition  of  time  (l.e.  £(x,t) 

•  EAn(x)exp(iont)) 

Velocities  in  the  x-,y-,  and  z-dlrections,  respectively 

Depth-averaged  velocities  in  the  x-  and  y-directions, 
respectively 

Salinity 

Turbulent  kinetic  energy 

Free  surface  elevation 

Interface  depth  of  two-layer  models 

Vertical  and  horizontal  turbulent  diffusion  coeffi¬ 
cients,  respectively 

Upper  and  lower  layers,  respectively 


*  From  Hamilton  and  Fucik  (1980). 
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Tablt  4 

Representative  Types  of  Wattr  Quality  Hodsls  of  Estusrlts* 


Reference 

Hydrodynamic  Circulation 
Submodel** 

Water  Quality 

State  Vsrlables 

Rarrett  and 
MoUovney  (1972) 

l-D  tidal ly  averaged 
box  model 

Organic  carbon,  organic 
nitrogen,  ammonia, 
nitrate,  DO 

Nlhoul  et  al. 

(1979) 

l-D  time-dependent 
dynamic  transport 
dispersion  model 

DO,  nitrate,  asmonla, 
iron,  manganese 

leendertsc  and 
Critton  (1971) 

2-D  (horizontal) 
time-dependent 
dynamic  transport- 
dispersion  model 

Conforms,  salinity  (chlo- 
rids*) ,  DO-ROD 

Chen  and  Or  lob 
(1972) 

.-D  time-dependent 
link-node  dynamic 
transport-dispersion 
model 

Teapsrature,  salinity, 
amonta  nltrogsn, 
nitrate  nitrogen, 
phosphorus,  DO- ROD, 
algae  (2),t  zooplank¬ 
ton,  fish  (2) 

DlToro  et  al. 

(1977) 

l-D  tldally  averaged 
transport-dlspsruton 
modsl 

DO-ROD,  phosphorus  (or¬ 
ganic  and  Inorganic), 
silica,  nlt/ogen  (ni¬ 
trate,  ammonia,  organic), 
phytoplankton,  zooplank¬ 
ton 

Kinetics  »»«*  Comments  Forcing 

Unnr  Tidal  amplitude, 

Ansroblc  kinetics  In-  rlverflow,  wind, 

eluded  when  DO  <  51  carbonacous  and 

saturation  nitrogenous  loads, 

temperature 

Thermodynamic  equilibrium  Tide,  rlvertlow 
assumed.  Redox  poten¬ 
tials  control  reaction 
rates  assuacd  to  be  due 
to  bacterial  activity 
Highly  polluted  -  no 
nutrient  limitation 
Linear 


Tide,  fresh  and 
stoiv  drain  flow, 
ROD  load,  wind 


Linear  and  Honod 


Tide,  rlvsrflow, 
waste  and  nutrient 
loads,  wind,  light 


Linear  and  Honod  Tide,  rlverflow, 

wind,  teuperature, 
light,  carbonecoua 
and  nitrogenous 
waste  loads 


Dahl-Kidsen 

(1978) 


box  models  or  1-0  time- 
dependent  transport' 
dispersion  model 


Selp  (1979) 
and  Selp  et  al. 
(1979) 


None  (observed  hydro¬ 
graphic  regime  used 
as  forcing  for  four- 
layer,  four  hortion- 
tal  section  model) 


Phytoplankton  (carbon,  nt-  Linear  and  Honod 
trogen,  phosphorus),  loo- 
plankton  (carbon),  detri¬ 
tus  (carbon,  nitrogen, 
phosphorus) ,  Inorganic 
nitrogen,  Inorganic 
phosphorus,  DO-ROD, 
sediment  carbon 


Rlverflow,  tide, 
temperature, 
light,  nutrient, 
and  organic  loada 


Benthic  algae  (8  age 

classes),  biomass,  zinc. 
Iron 


Model  of  population  Light,  salinity, 

dynamics  with  toxicity-  temperature,  nt- 

related  mortality  trogen  (limiting 

nutrient) 


*  Fro*  Hamilton  and  Fucik  (1980). 

**  l-D  •  one  c. senatorial;  2-D  •  two  dimensional, 
t  Numbers  In  parentheses  Indicate  the  number  of  simulations  available. 


Application 
Thamea  Estuary 
(U.K.) 


Scheldt 

(The  Netherlands; 


Jamaica  Ray 

(N.T.) 


San  Francisco  Ray 
Delta  System 


San  Francisco 
Ray-Delta  System 
Potomac  Estuary 


Various  Danish  fjords 
and  estuaries 
Prlnctpslly 
Ltmf Jorden 


Hardanger  fjord 
Sorf Jorden 
Trondhelms  [ jorden 
(Norway) 


Table  5 

iresentat Ive  Estuarine  Water  Quality  Models  Available  With  User  Documentation 
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In  accordance  with  letter  from  DAEN-RDC,  DAEN-ASI  dated 
22  July  1977,  Subject:  Facsimile  Catalog  Cards  for 
Laboratory  Technical  Publications,  a  facsimile  catalog 
card  in  Library  of  Congress  MARC  format  is  reproduced 
below. 


Hall,  Ross  W. 

Evaluation  of  marsh/estuarine  water  quality  and  ecological 
models  :  an  interim  guide  /  by  Ross  W.  Hall,  Jr. 
(Environmental  Laboratory,  U.S.  Army  Engineer  Waterways 
Experiment  Station).  --  Vicksburg,  Miss.  :  The  Station  ; 
Springfield,  Vn.  :  available  from  NTIS,  1981. 

18,  [6J  p.  ;  27  cm.  --  (Miscellaneous  paper  /  U.S, 

Army  Engineer  Waterways  Experiment  Station  ;  EL-82-1) 

Cover  title. 

"January  1982." 

Final  report. 

"Prepared  for  Office,  Chief  of  Engineers,  U.S.  Army." 

Bib! iogrnphy:  p.  12-18. 

1.  Estuaries.  2.  Marshes.  3.  Mathematical  models. 

4.  Water  quality.  1.  United  States.  Army.  Corps  of 
Engineers.  Office  of  the  Chief  of  Engineers.  II.  U.S. 

Army  Engineer  Waterways  Experiment  Station.  Environmental 
Laboratory.  Ill,  Title  IV.  Series:  Miscellaneous 
paper  (U.S.  Army  Engineer  Waterways  Experiment  Station)  ; 
EL-82-1. 
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